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Abstract 

Rail-guided vehicles (RGVs) are widely employed in automated freight handling system (AFHS) 
to transport surging air cargo. Energy-efficient routing of such vehicles is of great interest for 
both financial and environmental sustainability. Given a multi-capacity RGV working on a lin¬ 
ear track in AFHS, we consider its optimal routing under two-sided loading/unloading (TSLU) 
operations, in which energy consumption is minimized under conflict-avoidance and time window 
constraints. The energy consumption takes account of routing-dependent gross weight and dy¬ 
namics of the RGV, and the conflict-avoidance constraints ensure conflict-free transport service 
under TSLU operations. The problem is formulated as a mixed-integer linear program, and solved 
by incorporating valid inequalities that exploit structural properties of the problem. The static 
problem model and solution approach are then integrated with a rolling-horizon approach to solve 
the dynamic routing problem where air cargo enters and departs from the system dynamically in 
time. Simulation results suggest that the proposed strategy is able to route an RGV to transport 
air cargo with an energy cost that is considerably lower than one of the most common heuristic 
methods implemented in current practice. 

Keywords: Vehicle routing; energy efficiency; two-sided loading/unloading operations; pickup 
and delivery; conflict avoidance; time window constraints; multiple capacity 


1. Introduction 

Automated freight handling system (AFHS) is a type of automated material handling system 
(AMHS), which is widely adopted in facilities with massive material handling requests, such as 
freight terminals, distribution centers and production plants, to enhance system efficacy by mini¬ 
mizing operating cost and risk of human errors. As operating cost of AMHS can represent up to 
70% of the cost of a product HH2], it is critical to smartly design and operate AMHS to improve 
the overall economic and environmental performance. There has been a considerable growth of 
interest in studying such problems in both industrial and academic contexts. 
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This work considers improving AFHS installed in a freight terminal, which employs railed- 
guided vehicles (RGVs) to transport a tremendous amount of inbound and outbound cargo from 
their origins to destinations distributed along a linear track. The workload, which is especially 
high at its peak hour around midnight, leads to a conflict-prone environment that poses great 
challenges to terminal operations. The present AFHS has been developed to improve the termi¬ 
nal’s throughput while eliminating potential human errors. However, the design is not optimal 
especially in terms of energy efficiency. As energy consumption constitutes one of major sources 
of operating cost and has gained increasing attention for enabling a greener and sustainable earth, 
improving energy efficiency of AFHS and in particular developing an energy-efficient RGV routing 
strategy is of great interest to the industry. So far only heuristic methods have been employed 
to route RGVs in the current system, which motivates us to develop a rigorous mathematical 
programming method to improve the system performance. The new routing strategy also meets 
various service requirements such as delivery time windows (TWs) and avoidance of unloading 
deadlocks, etc., by incorporating them into the mathematical program from which the strategy 
results. This is contrast to existing methods which achieve that by abruptly compromising the 
system’s performance or by means of a posteriori sophisticated supervisory control |2j. 

1.1. Problem description 

A typical work area of the AFHS under consideration is depicted in Figure [lj An RGV is 
operated over a linear track to transport containers between work stations located along both sides 
of the track. Containers are queued at work stations, and will be picked up and delivered to their 
destined stations on either sides of the track by the RGV via so-called two-sided loading/unloading 
(TSLU) operations. The RGV and work stations are equipped with roller decks to support the 
TSLU operations. When the RGV is docked to a station, the roller decks rotate forward or 
backward accordingly to load or unload containers from or to either side of the track. The RGV 
can carry multiple containers subject to certain capacity limit. 

Each transportation of a container is initiated with a pickup and delivery (PD) request to 
the central control system. Midway drop-off is not allowed in the current practice because of the 
substantial overhead caused by frequent acceleration and deceleration of the RGV. Once a container 
is picked up, it remains on the RGV until being delivered to its destination. As containers enter and 
depart from the AFHS dynamically, the control system accumulates unfinished PD requests and 
aims at routing the RGV to pick up and deliver the associated containers in an optimal sequence 
so as to minimize energy consumption required for completing all transport requests subject to 
service quality and conflict-avoidance constraints for smooth operations. 
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Figure 1: Typical work area of an RGV in AFHS with exemplary PD requests. 

1.2. Related literature 

Sequencing PD tasks to be handled by a single vehicle is referred to as vehicle routing or job 
sequencing mm, and can be treated as a pickup and delivery problem (PDP) [Bjl?]. The problem 
is NP-hard in general due to a complex combinatorial nature, and it includes our problem under 
consideration as a sophisticated case. Our problem is further complicated by dynamic arrivals of 
PD requests. So far there have been a variety of research investigations on static/dynamic vehicle 
routing problems (VRPs) for different applications [8] |9l [TO]. The literature confined to a single 
vehicle is briefly reviewed as follows. 

Static VRPs assume that all transport requests are known a priori. Atallah and Kosaraju nu 
proved that the problem is polynomial-time solvable if the vehicle has unit capacity while confined 
to a linear track when no precedence and TW constraints are imposed on the transport requests. 
The same problem turns out to be NP-hard if the track topology changes to be a tree (or general 
graph) [121 ITO] . If the vehicle has multiple capacity, the problem is NP-hard even if the track is 
a simple path M- Another closely related problem, the PDP (for goods transportation) or the 
dial-a-ride problem (DARP, for passenger transportation) which includes TW constraints on PD 
requests has been well studied. Algorithms based on branch-and-cut or column generation are 
available for solving the problem of a small to medium size 0115]. Other related literature is on 
request/job sequencing in automated storage and retrieval systems (AS/RS), as referred to [1] for 
a recent review. 

The aforementioned literature all considered PDPs without loading constraints. In many ap¬ 
plications, however, constraints also appear on loading m- In the traveling salesman problem with 
pickup and delivery and LIFO loading (TSPPDL), a single vehicle must serve paired PD requests 
while both pickup and delivery must be performed in LIFO (last in first out) order. Heuristic 
algorithms for solving this problem were introduced in mm, and exact formulations and solu¬ 
tions were reported in [19, 20] which rely on tailored branch-and-cut algorithms. Another related 
problem is the traveling salesman problem with pickup and delivery and FIFO loading (TSPPDF) 
in which both pickup and delivery must be performed in FIFO (first in first out) order. Heuristic 
algorithms were introduced in [21] for solving the problem, and exact solutions were explored by 
using additive branch-and-bound [19] and tailored branch-and-cut [22|, respectively. 
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In practice, VRPs are dynamic in nature as PD requests appear stochastically in time. Research 
work on dynamic VRPs assumes the requests arrive in an unknown or known stochastic process 
Ei mg. With an unknown arrival process, deterministic algorithms have been proposed to treat 
dynamic VRPs and their performances are evaluated by competitive analysis (which compares the 
worst-case cost achieved by the algorithm for a dynamic problem and the optimal cost obtained 
for its counterpart having all requests being known a priori) |23 , Ml 15], However, the algorithms 
and their analyses all rely on the assumption that the dynamic VRPs do not have side constraints 
like loading, precedence and TW constraints. On the other hand, if the arrival process of PD 
requests is known, strategies based on sampling or Monte Carlo simulations are often used to 
handle dynamic VRPs [9]. Rolling-horizon based algorithms can also be employed if deterministic 
estimates of future requests are available |8]. Alternatively, heuristic algorithms may be developed 
by exploiting structural properties of a problem, as referred to [9] for a relevant review. 

In particular, among existing literature, e.g., [25U26U27U28], investigating operational aspects 
of AFHS (whereas investigations on higher-layer issues can be referred to |29j and the references 
therein), {28J studied a most relevant but simpler problem, in which the RGV is assumed to 
have unit capacity and the goal is to minimize RGV’s travel distance for completing all transport 
requests. The problem is a dynamic PDP with TW, FIFO queuing, and PD precedence constraints, 
and the investigation forms a pilot study towards obtaining a model and solution for the more 
complex routing problem considered in this work, of which partial results were reported in m- 

This work differs from existing literature in several aspects. Firstly, the problem under in¬ 
vestigation is a capacitated PDP under unique conflict-free service constraints. We completely 
characterize these service constraints under TSLU operations, which include the well-known LIFO 
and FIFO service constraints as two special cases. This is the first time that such kind of char¬ 
acterization has become available for transport service under TSLU operations, to the best of our 
knowledge. Secondly, the RGV routing problem aims to minimize total energy consumption of 
operating an RGV for completing all PD requests, which meets well with the interest of saving 
energy and reducing carbon emission in practice. This differs from the existing literature where 
travel distance or makespan (i.e., task completion time) is taken as the objective [6lE?, 31). Thirdly, 
structural properties of the new problem are exploited to reduce the problem domain and derive 
useful valid inequalities for improving the computational efficiency. Fourthly, a rolling-horizon 
approach is developed for treating the dynamic RGV routing problem, where a way of handling 
non-zero initial conditions (i.e., the RGV starts with nonempty load) is introduced and explained 
in detail. The new issues revealed and solved all explain the challenges of enabling energy-efficient 
routing of an RGV confronted in a real AFHS. 

The rest of the work is organized as follows. Section [2] perceives the static RGV routing problem 
as a sophisticated PDP problem, characterizes its conflict-free service requirements, and develops 
a full optimization model for it. Section [3] reformulates the initial model into a mixed-integer linear 
program (MILP), reduces its domain by removing infeasible solutions, and derives valid inequalities 
from the problem structure for expediting the solution process. Section [4] presents a rolling- 
horizon approach to treat the dynamic problem, as followed by comprehensive computational 


4 


studies performed on random instances in Section [5] Finally, conclusions are drawn in Section 
[6j Supporting materials that are helpful to understand the main results are collected in the 
Appendices. 


2. Static Routing Problem Formulation 


A multi-capacity RGV is routed to serve n PD requests with TW constraints on the deliveries. 
As midway drop-off is not allowed in current practice, fulfilling a PD request is equivalent to 
completing two tasks, pickup task (i.e., loading a container from its origin) and delivery task (i.e., 
unloading the container at its destination). Therefore, serving n PD requests is equivalent to 
completing 2 n tasks and the optimal routing solution is a processing order of the 2 n tasks for the 
RGV to consume least energy under service quality and feasibility constraints. 

To facilitate problem formulation, we assign a unique integer ID to each task. A pickup task 
is associated with an integer i, where 1 < i < n, and its corresponding delivery task associated 
with i + n. So the ith PD request means a pair of directed tasks i i + n, where ---> means a 
simple path that may contain multiple arcs. We denote the sets of pickup and delivery tasks as P 
and D , respectively, i.e., P = {1, 2,..., n} and D = {n + 1, n + 2,..., 2n}. In the example shown 
in Figure [lj there are seven PD requests indicated by dashed arrows, in which P = {1, 2,..., 7}, 
D = {8,9, ...,14} and the PD requests correspond to seven pairs of PD tasks {1 —-> 8,2 ---» 
9,..., 7 —■* 14}. For convenience, we refer to the PD requests by their pickup tasks, as P. 

The static RGV routing problem is then modeled by a directed graph G = (V, A), where a 
vertex i in V represents a pickup or delivery task with the ID of i and an arc (i, j) in A represents a 
plausible processing order between the two tasks i and j, i.e., whether the task j could be processed 
right after the task i. Two virtual tasks, 0 and 2 n + 1, are added to represent the start and end 
positions of the RGV, respectively. Define V' = PUD = {1,2,..., 2n} and V = {0}UV / U{2n+1}. 
The arcs (i, j) for all i, j G V' and i,j G V give rise to arc sets A' and A, respectively (both of 


which can be reduced by removing infeasible arcs in Section 3.2). 


A feasible routing solution is then a path starting from vertex 0, going through all vertices 
in V' exactly once and ending at vertex 2n + 1. While a path can generally be formulated by a 
sequence of binary decision variables Xij which indicate the arcs (i,j) used in the path, another 
group of binary decision variables yfj, which indicate the arcs traversed during the service of each 
request k is also required for enforcing conflict-free service constraints as revealed next. The new 
problem and its unique model differentiate our work from existing literature. 


2.1. Conflict-free service constraints 

In practice there are physical constraints, such as containers on an RGV cannot swap positions, 
that restrict the TSLU operations. The multi-capacity RGV thus must be routed to load and 
unload containers in a valid order in order to avoid infeasible operations and deadlock. Here, 
infeasible operation means an operation that attempts to unload a container before another one 
which is however impossible to realize in practice for physical constraints; and deadlock means that 
a container cannot be unloaded without temporally dropping off another container and it remains 
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so if the order of unloading the two containers is reversed. The requirements are equivalent 
to imposing conflict-free service constraints on the TSLU operations. Before describing these 
constraints, we classify all PD requests to be handled by the TSLU operations into four types: 

P 1 = {i e P : ch = ai+ n = 0}, P 2 = {i E P : <h = a i+n = 1}, 

P 3 = {i E P : a,i = 0, a l+n = 1}, P 4 = {i E P : a % = 1, ai +n = 0}, (1) 

where a* is an indicator variable to show on which side task i occurs, and it indicates the north side 
if at = 0 and the south side otherwise. The four types of PD requests are illustrated in Figure [2| 

For the example shown in Figure [lj we have P 1 = {1}, P 2 = {3, 5} P 3 = {2, 6} and P 4 = {4, 7}. 

Type-1 Type-2 Type-3 Type-4 

Figure 2: Four types of PD requests. 

Various conflicts may arise if tasks belonging to different types of PD requests are processed 
sequentially. The types of conflicts reviewed in the literature (e.g., [201 [221 HI]) can be regarded 
as loading/unloading systems restricted to one type of the PD requests defined above. Thus, their 
conflict-free service constraints can be satisfied by sticking to a certain task service order: If PD 
requests are all of Type-1 (or Type-2), the LIFO service order will suffice; and if PD requests are 
all of Type-3 (or Type-4), the FIFO service order will do. 

Under TSLU operations, however, a single LIFO or FIFO service order cannot guarantee 
conflict-free services. As there are four mixed types of PD requests, we will face 10 groups of 
scenarios for serving different combinations of PD requests. Among them, one group is always 
conflict-free in which tasks of different types of requests can be served in any order, and that the 
rest 9 groups have to be constrained to avoid conflicts. Specifically, the conflict-free scenarios 
are concerned with service of PD request pairs P 1 ~ P 2 , and the rest 9 groups of scenarios are 
concerned with service of PD request pairs P 1 ~ P 4 /P 3 /P 4 , P 2 ~ P 2 /P 3 /P 4 , P 3 ~ P 3 /P 4 , 
p 4 ~ p 4 ; where the symbol / means “or”. Different constraints may be imposed in the 9 groups 
of scenarios to ensure conflict-free services. 

The ten groups of scenarios and their associated service constraints can be classified into six 
cases below. The hrst two cases correspond to conventional scenarios with LIFO or FIFO loading 
restrictions, and the next three cases are unique to the routing problem under consideration, and 
the last case describes conflict-free scenarios where no special service constraint is required. 

2.1.1. Case 1. P 1 ~ P 1 and P 2 ~ P 2 : LIFO service 

As illustrated in Figure [3j LIFO service order should be maintained between any pair of PD 
requests in P 1 (or P 2 alike) for conflict avoidance. Given requests j. k E P 1 which share the RGV 
for a certain period of time, if the pickup task j is processed after the pickup task k, then the 
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delivery task j + n must be processed before the delivery task k + n. 




Figure 3: LIFO service: P 1 ~ P 1 and P 2 ~ P 2 . 


To characterize this mathematically, we introduce binary decision variables y\ 3 for all request 
k E P and arc (i. j) in a feasible set. We have y\ 3 = 1, if arc (i,j) is traversed on the path from 
vertex k to vertex k + n (i.e., during the service for request k); and y\ 3 = 0 otherwise. The LIFO 
service order is enforced by the following constraints: 


£ y% 


£ 


v k - 


i:(i,j+n)(E:A' 


Vj E € P\ 

Vj E P 2 \{k},k€ P 2 . 


( 2 ) 


Constraint means that if task j is processed between k and k + n, then task j + n must also 
be processed between k and k + n, namely, the LIFO service order is implemented. 


2.1.2. Case 2. P 3 ~ P 3 and P 4 ~ P 4 ; FIFO service 

Similarly, FIFO service order should be maintained between any pair of PD requests in P 3 
(or P 4 alike) for conflict avoidance, as illustrated in Figure [3j This requirement can be met by 
enforcing constrain ([ 3 ]) below, which means that if a pickup task of request k in P 3 (or P 4 ) is 
processed before the service of request j of the same type, then the delivery task of request k must 
be processed before the completion of request j, i.e., the FIFO service order is enforced. 



X/ y ij + yh+n < 

i:(i,j)EA' i:(i,j+n)(E:A' 


\/j E P 3 \{k},k E P 3 , 
Vj E P 4 \{k},k € P 4 . 


(3) 


2.1.3. Case 3. P 1 ~ P 3 and P 2 ~ P 4 : crossing first in (CFI) service 

This case is illustrated in Figure [ 5 ] Assume j E P 3 , k E P 1 and they are simultaneously served 
by an RGV in a certain time interval. The PD request j is a crossing request whose pickup and 
delivery locations are on opposite sides of the track. If the pickup task j is processed after the 
pickup task k but before the delivery task k + n, then the delivery task j + n cannot be performed 
because the container associated with the request k will block the way out. Thus, in this case, the 
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crossing request must be handled first for ensuring conflict-free service, leading to the so-called 
CFI service order. This applies to the case for j E P , k E P 2 alike. 



Figure 5: CFI service: P 1 r \j P 3 and P 2 


rsj 


P 4 . 


The CFI service order can be enforced by the following constraints: 


Vij = 0E A', 


Vj E P 3 , k E P 1 , 
Vj E P 4 , A: E P 2 . 


(4) 


Also note that, for j E P 1 , k E P 3 (or j E P 2 ,A; E P 4 ), we can obtain another set of conflict-free 
service constraints, Yhi:(i,j+n)eA'Vi,j+n — Y!v.(i,j)eA' Viji meaning that if task j + n is processed 
between tasks k and k + n, then task j must also have been processed. As shown in Appendix [Aj 
this set of constraints are equivalent to the constraints in Q and hence can be ignored. 

2.1.4- Case 4 • P 1 ~ P 4 and, P 2 ~ P 3 / crossing last out (CLO) service 

Similarly, CLO service order should be maintained in Case 4 as illustrated in Figure [6j Assume 
that j E P , A: E P 1 and they are simultaneously served by the RGV in a certain time interval. 
Request j is still the crossing request. This time, the pickup tasks j and k can be processed in a 
flexible order, but the delivery task j + n must be processed after the delivery task k + n because 
the container associated with request j is queued behind the one associated with request k on the 
RGV. 



The CLO service order can be enforced by the following constraints: 


Di,j+n ~ j + n ) ^ At i 


Mj e P 4 ,k e P 1 , 
Vj e P 3 ,fc e P 2 . 


(5) 


As in the CFI case, for all j E P l ,k E P 4 (or j E P 2 ,k E P 3 ) we can obtain another set of 
conflict-free service constraints, ^2i-.u,j)eA' Vij - Y^i:(i,j+ n )<=A' Vi.j+m which are equivalent to the 
constraints in ([5j and hence can be ignored. A proof of this fact is again referred to Appendix [Aj 


2.1.5. Case 5. P 3 ~ P 4 : Deadlock 


This case corresponds to the deadlock case for any pair of PD requests P 3 ~ P 4 as illustrated 
in Figure [7j Such pair of PD requests should never be simultaneously served by the RGV because 
the corresponding containers will block their way out from each other. 



Figure 7: Deadlock: P 3 ~ P 4 . 

The deadlock can be avoided by enforcing the following constraint, which completely avoids 
overlap between services of the aforementioned two types of PD requests: 



It can be shown that the above condition implies that yff ' - +n = 0, for all ( i,j + n) E A' with j, k in 
the domains given in ([6]). 

2.1.6. Case 6. P 1 ~ P 2 : Free case 

This case corresponds to the free case illustrated in Figure [8j in which any pair of PD requests 
do not interfere from each other. The requests can thus be freely served by the RGV with no 
dedicated service constraints. 



Figure 8: Free service: P 1 ~ P 2 . 


Remark 1. The LIFO and FIFO service constraints formulated in 0 and 0 have forms con- 
ciser than those developed in m where two additional groups of binary variables with the same 
dimensionality of y - are required. The current formulation avoids those binary variables and the 


constraints associated, and hence is computationally more efficient in general. 

2.2. Problem formidation 

This subsection develops a mathematical model of the energy-efficient RGV routing problem. 
The major parameters and notations used are listed below for ease of reference: 
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(a) Container with unit load. (b) Containers with two units of load. 
Figure 9: Containers with different capacities. 


n 

P 

D 

V 
A' 

V 
A 
Q 
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Si 

[ e O h] 
nj 

tij 


The number of PD requests; |P| = \D\ = n; 

The set of pickup tasks, P = {1, 2,..., n}; 

The set of delivery tasks, D = {1 + n, 2 + n,..., 2n}; 

The union of pickup and delivery tasks, V' = P U D ; 

The set of arcs (i, j) with i,j £ V'• 

The full tasks including virtual start and end tasks, V = {0} Ui"U {2n + 1}; 

The set of arcs (i,j) with i,j £ V; 

The RGV’s capacity in units of load (and every unit has the same weight j^J); 

The units of load associated with task i, satisfying qi > 0, qi+ n < 0 and qi = —qi+ n 
for all i £ P. There are several types of containers with different sizes. For example, 
the container in Figure [9 (a) | admits one unit of load, and so qi = 1, qi+ n = —1; and the 
container in Figure 9(b) admits two units of load, and so qi = 2 ,qi+ n = —2; 

The operating duration of task i £ V\{2n + 1}, with so = 0; 

The delivery TW associated with each PD request i £ P; 

The RGV travelling distance between tasks i and j; 

The RGV travelling time between tasks i and j\ 


There are four groups of decision variables: 


Xij The binary variable to indicate whether task i is processed right before task j; 

Dij The binary variable to indicate whether the arc (i,j) is traversed on the path from 
vertices k to k + n, where k £ P; 
bi The time starting to handle task i £ V, with b o = 0; 

Wi The units of load on the RGV upon completion of task i £ V, with wo — 0. 

The static optimal RGV routing problem can be formulated as a mixed integer program defined 
in 0-©- For convenience of presentation, we term it as the pickup and delivery problem with 
TSLU operations, or PDP-TSLU for short. 


1 Generalization can be made to associate each unit of load with a different weight if such information is available. 
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PDP-TSLU : 

min ^2 Cij(wi)xij 

(7) 




subject to, Xij = 1 

Vi G V\{2n + 1}; 

(8) 

j: ( i,j)eA 



J2 x i i = 1 

Vi € V\{0}; 

(9) 

j: (, j,i)£A 



bi-\-n A bi "f- Si 4“ t-i^i-i-n 

Vi G P; 

(10) 

bi®l > bi + Si 

V&i©i /0,RP; 

(11) 

J bj > bi + Si + tij 
lj | Wj >Wi + qj 

V(i, j) G A, j G V’ U {2n + 1}, 

V(i, j) € A, j £ V'; 

(12) 

e% F bi-\. n Si-\- n ^ li 

Vi G P; 

(13) 

max{0, q,} < Wi < min{Q, Q + gj 

Vi G V'; 

(14) 

(1 if * = k 



— ijji = < — 1 if * = k + n 

Vi G V', k G P; 

(15) 

j:(j,i)eA' otherwise 



©-©; 


(16) 

Vij < Xij 

V(i, j) G L, 1 G P; 

(17) 

Vij e {0, 1} 

V(i,i)Gi',*GP; 

(18) 

Xij € {0, 1} 

V(i, j) G A; 

(19) 

bi > 0, Wj > 0 

ViG V' U{2n + 1}, j G V'. 

(20) 


The objective function in ([7]) measures the total energy consumption and is a bilinear function 
in Wi and x t j , whose specific form is derived as follows. Each arc (i,j) G A is associated with 
an energy cost Cij(wt ) (as explicitly depends on the load weight) and a travel time Uj. With the 
operating profile of the RGV shown in Figure |T7)) the arc travel time can be expressed as an explicit 
function of the arc distance rij as 


Uj — 



if 0 < Tjj < 2ri, 
if rij > 2n, 


( 21 ) 


where the parameters a, v c , t\ and r\ are defined in Figure [10J Let wrgv be the net weight of 
the empty RGV and Wi denote the load weight upon leaving vertex i. The energy cost Cij(wi ) is 
calculated as the work done by the RGV for traversing the arc (i, j). By Newton’s law we have, 


c ij ( w i) — 


Hgrij(wRGV + Wi) 

< arijiwRGv + Wi) 

^ (2(a - fig)n + fj.grjj) (■ w RG v + Wi) 


if a < gg- 

if a > fig and 0 < r^ < 2ri, 
if a > fig and > 2r\, 


( 22 ) 


where fj is the rolling friction coefficient while the RGV moves on the track and g is the gravitational 
acceleration, both of which are given constants. Although fj , g and r ? y are given constants, the 
energy cost Cij{wi) are variables controllable by adjusting the task service order. This cost is 
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• RGV Specifications 

- Traveling velocity (v c ) = 3 m/s 

- Acceleration/deceleration (a) = 0.6 m/s 2 

• Standard Time-Graph Calculations 

- v = at 

- r = 0.5 at 1 


• To compute acceleration/deceleration 
time/distance 

-1 = via 
-r = 0.5al 

• Specifically, 

-h = h ~ h = 5 s 
-r 1 = r 3 -r 2 = 7.5 m 


Speed ^ 
3 m/s 


0 1 1 


Speed ^ 
3 m/s 


t 2 t 3 Time (s) 


0 n 


r 2 r 3 Distance (m) 


Figure 10: Operating profile of an RGV. 


generally increasing in Wi, in contrast to conventional routing cost which is merely determined by 
the arc distance . In the meanwhile, we set Cj^n+i = 0 for all i € D, to account for the practice 
that the RGV stops at the last delivery. 

We briefly explain the constraints of PDP-TSLU. Constraints ([8]) and Q ensure that each 
task is served exactly once, with the service starting and ending at virtual tasks 0 and 2n + 1, 
respectively. Constraint © states that each pickup task should be processed before its paired 
delivery task, and constraint © means that containers at the same work station are picked up 
in FIFO order. In the term i ® 1 represents the pickup task queued right behind task i at 
the same station. For the example shown in Figure [lj we have 1 © 1 = 2 and 4 © 1 = 5 . 

The indicator constraint (12) ensures the consistency of service time and load variables and 
also the route connectivity, which can be linearized using the big -M method (refer to Appendix 
[B]for details). In (12), the load inequality constraints are valid alternatives of the corresponding 
equality constraints for the optimization considered, but this is true for the time constraints if 
only latest but earliest service constraints are active. Constraint (13) imposes a TW constraint 


on completion time of each request. Constraint (14) enforces the RGV’s capacity limit, which is 
often tighter than applying the obvious bound 0 < w 2 < Q for all i £ V'. 

Constraint (15) ensures a path from vertices k to k + n for any PD request k that does 


not pass through the virtual start or end vertex. Constraint (16) consists of unique conflict-free 
service constraints resulting from the TSLU operations as revealed in the previous subsection, and 
constraint © binds the overall routing decisions with the individual routing decisions of every 
PD request. Finally, constraints (18)-(20) specify domains of the decision variables. 


Remark 2. If the objective function is replaced with total travel distance or time, the PDP-TSLU 
reduces to TSPPDL confined on a linear track if there are only PD requests of Type-1 (or -2), and 
to TSPPDF confined on a linear track if there are only PD requests of Type-3 (or -f). In the two 
reduced cases, the conflict-free service constraints 0-0 collapse to mean the LIFO and the FIFO 
service order, respectively. 
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3. Solution approach 

This section reformulates the bilinear PDP-TSLU into an MILP, and then reduces the arc set 
(and so the variables and constraints associated) based on delicate insights into the PDP-TSLU, 
and next introduces a couple of valid inequalities that exploit structural properties of the problem 
for solving it more efficiently. 

3.1. Reformulating PDP-TSLU as an MILP 

Let aij = arij and fii 3 = 2(a — g.g)r\ + ^grij. The PDP-TSLU is bilinear in its present form 
because the objective function has the following explicit form: 

Y Cij(wi)xij = Y a U (wrgv + Wi)xij + Y Pij (wrgv + Wi) Xij, ( 23 ) 

( i,j)eA (i,j)eA,r i: j<2n (i,j)eA,nj>2n 

which contains bilinear terms of WiXij. This makes the problem difficult to optimize directly. 
Fortunately, the objective function can be linearized by introducing new decision variables and 
additional linear constraints. 

Without loss of generality, let us assume a > fj.g (which means that braking is required for 
stopping a moving RGV at a deceleration of a). In this case, the objective function of the PDP- 
TSLU is equivalent to 


Y c ij( w i) x ij = Y Zi ' 

(i,j)(zA ieV\{2n+l} 

subject to, for all feasible index j, 

, | aij (wrgv + Wi) V(i, j) e A Tij < 2 n, 

Xij = 1 =>• Zi = < 

[Pij(w RGV + V(z, j) e A, rij > 2n. 

Note that it is unnecessary to introduce Zij for each x VJ , because there always exists one and only 
one immediate successor j to i such that = 1 by constraint ([8]). 

The above indicator constraints can be handled directly by solvers embedded in software like 
CPLEX [32). Alternatively, they can be reformulated into linear forms as 

Zi > aij (wrgv + Wi) - 7 ij(l - Xij), V(i,j) € A, nj < 2 n, 

Zi > Pij(w RGV + w^ - 7^(1 - x^), V(i, j) e A, > 2n, 

where 7 ij are upper bounds of aij(w RG v + Wi) if nj < 2 ri and of fiij{w RG y + Wi ) if nj > 2 ri- 
It is feasible to set 7 ^ = a t j(w RG v + min{Q, Q + %}) and 7 ^ = /% ( W RC!V + rriin{Q, Q + ^}) 
for the two cases, respectively. Note that the inequality relaxations of the indicator constraints 
do not change the optimal solution because the objective function minimizes the sum of z t wtih 
unit coefficients. Consequently, the PDP-TSLU becomes an MILP solvable by standard solvers. 
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The solution process can further be enhanced by exploiting structural properties underlying the 
problem, which are pursued in the next two subsections. 


3.2. Arc set reduction 

For practical reasons, not all arcs of the complete graph are feasible. By exploiting structural 


properties of the problem, the arc set A in PDP-TSLU is reducible to the one defined in (24), 


where the relation i<\j (or i O j ) means that pickup task i is queued in front of (or behind) pickup 
task j at the same station and the term i © 1 was defined when explaining ©■ 


A = j(0, j) : j G P such that {k : k <] j} = 0 j 

U j(z, 2n + 1) : i G D such that {k : k\> i — n} = 0 j 
{(*, j) '■ i £ P, j £ P k) D\ U? =1 Sp{\ Uf =1 S l D l | 
{(*, j) '■ i £ D, j E P (J D\ uf =5 S’ Pl \ uf =4 S l Dl | 
i, j) : i,j G P U D such that i = j or i = j + n j, 


U 


U 


(24) 


where 


S i Pjl ±{j€P:j<i}, 

Sp, 2 = {j £ P '■ jA> i © 1 } , 

Sp, 3 — {j £ P '■ j £ P 3 i if * 6 P 1 U -P 4 } , 

{j e p ■. j e p\iii e P 2 u P 3 }, 

Sp 5 - i:i e p '■ 3 

Sp^ 6 = {j G P : {k G P 3 : i — n <1 k < j} ^ 0, if i — n G P 1 } 

Spj = {j G P : {k G P 4 : * — n <1 k <1 j} ^ 0, if z — n G P 2 }, 

8 — {j € P : | {& € -P 3 : z — n < k <1 j} | > Q, if i — n G P 3 } 

Pp ,9 — {j G P : | {& G P 4 : z — n < k <3 j} | > Q, if i — n G P 4 }, 

5b, l - {j e D ■ 

5^2 4 {j g D : j - n G P 1 U P 4 \{z}, if i G P 1 U P 3 } 

5b ,3 = {j G P : j - n G P 2 U P 3 \{z},if z G P 2 U P 4 } , 

Sb ,4 — {j GD: j — n > i — n, if z — n G P 1 U P 2 } , 

5b, 5 — {j GO: j — n <\ i — n, if z — n G P 3 U P 4 } , 

Pb,6 — {i £ P : j — n £ Sp 8 U Spg} , 

5b i7 = {j G D : j - n G P 2 U P 4 , if z - rz G P 3 } , 

5b, 8 = {j G D : j - n G P 1 U P 3 , if z - n G P 4 } . 

Of the reduced arc set A, the first subset consists of arcs connecting the start vertex with the 
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head vertex of each pickup queue. The second subset consists of arcs connecting each delivery 
vertex, whose paired pickup vertex is the tail of a pickup queue, with the end vertex. The third 
subset consists of arcs connecting a pickup vertex with all other pickup vertices that are neither 
predecessors of the current vertex (S p x ) nor successors other than the nearest one in the same 
pickup queue (S P 2 ) and meanwhile services of which avoid unloading deadlock (i.e., excluding 
Sp 3 U S P4 ), and with delivery vertices whose paired pickup vertices are not successors of the 
current vertex (S l D x ) and meanwhile services of which avoid unloading deadlock (i.e., excluding 

^b, 2 u ^b, 3)- 

The fourth subset of A is most complicated. It includes arcs connecting a delivery vertex with 
pickup vertices which are neither predecessors of the pickup vertex currently in pair (S p 5 ) nor 
successors that queue after the first successor of which the service leads to unloading deadlock 
(S p 6 U S p 7 ), nor successors that queue after the pickup vertex currently in pair for more than 
Q request distance away while the requests are of Type-3 or Type-4 as the same as the current 
request (S p 8 U S l p 9 ), and also includes arcs connecting a delivery vertex with delivery vertices 
whose paired pickup vertices are neither successors of the pickup vertex currently in pair if the 
current request is of Type-1 or Type-2 (S l D 4 ) nor predecessors or successors that queue after the 
pickup vertex currently in pair for more than Q request distance away if the current request is of 
Type-3 or Type-4 (S l D 5 U S l D 6 ), and with delivery vertices that are not on the same side of the 
current vertex while the current request is of Type-3 or Type-4 (i.e., excluding S l D 7 U S l D 8 ). 

Overall the above reductions of the arc set owe to the precedence, loading and capacity con¬ 
straints inherent in the transport service. To have a feel of the power of the above reductions, 
consider the toy example shown in Figure [lj The arc set is reduced from an original size of 225 
(15 2 , for the complete graph) to a size of 114 (for a much sparser graph), cutting almost half of 
the full arcs. A direct consequence of such reductions is that decision variables and constraints 
associated with the identified infeasible arcs are eliminated before formulating and solving the 
problem, contributing to a conciser model and hence higher computational efficiency in general. 
It should be noted, however, that the elimination of infeasible arcs is incomplete, which bears the 
necessity of using the precedence, loading and capacity constraints in PDP-TSLU to achieve that. 

3.3. Valid inequalities for PDP-TSLU 

For a branch-and-cut algorithm to solve the PDP-TSLU, valid inequalities exploiting struc¬ 
tural properties of the problem can be useful to tighten the lower bounds and avoid unnecessary 
branching which will consequently expedite the process of searching for an optimal solution. 

The PDP-TSLU is essentially a generic PDP subject to new loading/unloading constraints. 
Therefore some valid inequalities known for a generic PDP apply to the PDP-TSLU. New valid 
inequalities can also be derived from the unique conflict-free service constraints in a way similar to 
those from LIFO or FIFO loading constraints as presented in |201 [22] . Although there are dozens 
of such valid inequalities, not all of them achieve good trade-offs between performance benefit 
and implementation cost: According to the numerical studies reported in [20L122]. some kinds of 
valid inequalities require complicated implementations (with respect to the demand for computer 
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memory and computational burden) but only slightly accelerate the solution process. Similar kinds 
of valid inequalities are therefore ignored for the PDP-TSLU without further validation. 

3.3.1. Valid inequalities inherited from generic PDP 

Let S be a subset of P U D. Define x(S) = Yhijes x iP an< ^ ^ 7r ('^) = {z £ P|* + n E S} and 
cr(S) = {i + n E D\i E S n P} which denote the sets of predecessors and successors of certain 
vertices in S, respectively, where i + n is the delivery task paired with the pickup task i E P. 
Define S = V\S. Then the following subtour elimination constraints are well-known for the PDP 
(as a generic property of TSP problems) |2(>, [22] ; 

x(S) < |5| - 1, V|S| > 2. (25) 

By taking account of precedence restrictions in PDP, these constraints can further be strengthened 
as mmm- 


(s) +22 22 x p 

+ £ 

22 Xij < S’ - 1, 

(26) 

i£S j€§nir(S) 

i£SC\7v(S) j£S\n(S) 


is) + 22 22 x p 

+ £ 

22 X V ^ I 5 ! - L 

(27) 

iE§r\a(S) 

ieS\a(S) jeSna(S) 



Another set of valid inequalities for the PDP are known as lifted and Df inequalities. 
The D £ and Df. inequalities were hrstly introduced by Grotschel and Padberg [33] and later 
strengthened by Cordeau [35]. The lifted inequalities apply to any ordered set S = ■ ■ ■, ik} V 


V for k > 3, and take the form of 



k k— 1 

k-1 h-1 


'22 Xi hih+1 + 2 22 x ihh + 

2222/ Xihii 22 x j+n,h < k — i, 

(28) 

h=1 h=2 

h= 3 1=2 j-\-n£Sr\a(S) 


k k—1 

k h—1 


22 Xi h j -h+i + ^ 22 x 

hih + 22i22i Xihii ^ _ i) 

(29) 

h=1 h =3 

h=A 1=3 j£Sr\ir(S) 



where ik+i = i\- 

3.3.2. Valid inequalities derived from LIFO and CLO service restrictions 

We derive new valid inequalities from the LIFO and CLO service constraints as unique to 
PDP-TSLU. Given i E P 1 , j E P 1 \{z}, if Xij = 1 is a feasible integer solution, then the pickup 
and delivery sequence must satisfy 0~<i~<j~<j + n~<i + n~<2n- (-1, where the operator -< 
represents the precedence relationship. To enforce LIFO service order under TSLU operations, 
the predecessor of i + n can only be a pickup task in P 2 U P 4 or a delivery task of a request in 
P\(P 4 U {*}) (including j). Meanwhile, a valid successor of j + n can only be a pickup task in 
P\{P 3 U {i,j}} or a delivery task of a request in P 2 U P 3 U {i}. Similar reasoning can be applied 
to determine the possible predecessors of i + n and successors of j + n for i belonging to P 4 and 
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for following CLO service order under TSLU operations. In summary, if arc ( i,j ) is included in 
the routing path, then the set of possible predecessors of i + n are given as 


V i+n (i,j) = P 2 U P 4 U {k + n : k G P\(P 4 U {*})}, Vi € P\j G P\{i}, 
V i+n {i,j) = ( P 2 U P 4 \{i}) U {k + n : k G P\(P 3 U {*})}, Vi G P 4 ,j G P 1 , 

iP i+n (i, j) = P 1 U P 3 U {A: + n : k G P\(P 3 U {i})}, Vi G P 2 , j G P 2 \{i}, 

P t+ n(i, j) = ( P 1 U P 3 \{i}) U {£: + n : k G P\(P 4 U {i})}, Vi G P 3 ,j G P 2 , 

and meanwhile the set of possible successors of j + n are 

S j+n (i,j)= (P\(P 3 U{i,j}))u{fc + n: k G P 2 U P 3 U {i}}, Vi G P 4 ,j G P 4 \{i}, 
Sj+n(i,j) = (P\(P 3 U {i, J})) U {/c + n : k G P\(P 3 U {j})}, Vi G P\j G P 1 , 

5,+n(i, j) = (P\(P 4 u {i, j})) U {fc + n : A; G P 1 U P 4 U {i}}, Vi G P\j G P 2 \{i}, 
5 i+n (i,j) = (P\(P 4 U{i,j}))u{fc + n: k G P\(P 4 U {j})}, Vi G P 3 ,j G P 2 . 

On the other hand, given i G P 1 , j G P 1 \{i}, if Xj + n,i+n = 1 is a feasible solution, then the 
pickup and delivery sequence again must follow the aforementioned precedence order. To enforce 
the LIFO service order under TSLU operations, the predecessor of j can only be a pickup task in 
{i}UP 2 UP 4 or a delivery task of a request in P\(P 4 U{i, j}). Meanwhile, a valid successor of i can 
only be a pickup task in P\(P 3 U {i}) or a delivery task of a request in P’llP 3 . Similar reasoning 
can be applied to determine the possible predecessors of j and successors of i for i belonging to P 4 
and for satisfying the CLO service order under TSLU operations. In summary, if arc (j + n,i + n ) 
is included in the routing path, then the set of possible predecessors of j are given as 

Vj{j + n, i + n) = {i} U P 2 U P 4 U {k + n : k G P\(P 4 U {i, j})}, Vi G P 1 , j G P 4 \{i}, 

Vj(j + n,i + n) = {0}UP 2 UP 4 U{fc + n: P\(P 3 U {i,j})}, Vi G P 4 , j G P 1 , 

Vj(j + n, i + n) = {i} U P 1 U P 3 U {k + n : k G P\(P 3 U {i, j})}, Vi G P 2 , j G P 2 \{i}, 

Pj (j + n, i + n) = {0} UP 1 UP 3 U{i , + n: P\(P 4 U {i, j})}, Vi G P 3 , j G P 2 , 

and meanwhile the set of possible successors of i are 

S t {j + n, i + n) = (P\(P 3 U {i})) U {fc + n : k £ P 2 U P 3 }, Vi G P 1 , j G P 1 \{i}, 

«Si(j + n, i + n) = (P\(P 3 U {i})) U {/c + n : k G P 1 U (P 4 \{i})}, Vi G P 4 , j G P 1 , 

5j(j + n, i + n) = (P\(P 4 U {i})) U {fc + n : /c G P 1 U P 4 }, Vi G P 2 , j G P 2 \{i}, 

<Sj(j + n, i + n) = (P\(P 4 U{i})) U {k + n : fc G P 2 U (P 3 \{i})}, Vi G P 3 ,j G P 2 . 

In the second and fourth cases of predecessors set of j , the vertex 0 instead of i is included in the 
set because j + n -< i + n does not imply i ~< j in those two cases. 

With the above insights, four groups of valid inequalities can be obtained for the PDP-TSLU, 
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each of which indicates a group of arcs that are incompatible with the arc (i. j ) or (j + n,i + n ) as 
included in a valid routing path. 


Proposition 1. For each of the vertex pairs ieP 1 U P 4 , j £ P 1 \{i}, or i £ P 2 U P 3 , j £ P 2 \{z}, 
t/ie following inequalities hold for the PDP-TSLU: 




^ ( Xlj+n ^ I; 

(30) 



''f ] Xj+n,l — 1) 

(31) 




(32) 

l$Vj(j+n,i+ri) 


y xu < i, 

(33) 

l£Si(j+n,i+n) 



where the variables are zero if their corresponding arcs are not defined in A. 

The next inequality comes from the fact that whenever an arc (i, i + n) for any i £ P 3 U P 4 
is used, it follows from the CLO loading restriction that the vehicle must arrive empty at i and 
leaves empty from i + n. 

Proposition 2. For each of the vertex pairs i £ P 4 , j £ P 1 , or i £ P 3 , j £ P 2 , the following 
inequality holds for the PDP-TSLU: 


Xij T Xji + 3'j,j+n T j+n 1; 

where the variables are zero if their corresponding arcs are not defined in A. 


(34) 


The following inequality merges the two valid inequalities + Xi+ n j+ n + Xj+ n ,i < 1 and 
x^ + Xi +n j +n + Xi +n j < 1, for each pair of i and j in proper sets. 

Proposition 3. For each of the vertex pairs i £ P 1 U P 4 , j £ P 4 \{i}, or i £ P 2 U P 3 , j £ P 2 \{i}, 
the following inequality holds for the PDP-TSLU: 


Xij T Xi-\. n j+n T Xj-\. n i -\- Xi-\-n,j T 1, 

where the variables are zero if their corresponding arcs are not defined in A. 


(35) 


3.3.3. Valid inequalities derived from FIFO and CLO service restrictions 

This subsection derives new valid inequalities pertinent to the FIFO and CLO service con¬ 
straints of the PDP-TSLU. Given i £ P , j £ P 4 \{?'}, if x^ = 1 is a feasible solution, then the 
pickup and delivery sequence must satisfy 0-<i~<j~<i + n~<j + n-<2n + l. To enforce the 
FIFO service order under TSLU operations, the successor of i + n can only be a pickup task in 
P\(P 3 U {?', j}) or a delivery task of a request in P 2 U {j}. Meanwhile, a valid predecessor of j + n 
can only be a pickup task in P 2 U (P 4 \{i,j}) or a delivery task of a request in P 1 U P 2 U {i}. 
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Similar reasoning can be applied to determine the possible successors of i + n and predecessors 
of j + n for z belonging to P 1 and for satisfying CLO service order under TSLU operations. In 
summary, if arc (i,j) is included in the routing path, then the set of possible successors of i + n 
are obtained as 


S i+n (i,j) = ( P\{P 3 U {i,j})) U {k + n: k E P\(P 3 U {z})} , 
S i+n (i,j) = (.P\(P 3 U U {k + n: k E P 2 U {j}} , 

S i+n (i,j) = (■ P\(P 4 U U{k + n:ke P\(P 4 U { z })} , 

Si+n(i,j) = {P\(P 4 U{z, j})) U {k + n : k E P 1 U {j}} , 


VzEP\j EP 4 , 
VieP\j eP 4 \{i}, 
VieP 2 ,jeP :3 , 
VzEP 3 ,jEP 3 \{z}, 


and meanwhile the set of possible predecessors of j + n are 


Pj+n(iJ) = P 2 U (P 4 \{j}) u {fc + n : fc E P\(P 3 U {j})} , 
V j+n (i,j) = P 2 U(P\{i,j})u{k + n : kP'uF'uji}}, 
^■+n(b j) = P 1 U (P 3 \{j}) U {/c + n : fee P\(P 4 U {j})} , 
Vj+n(i,j) = P 1 U (P 3 \{i, j}) U {fc + : k G P 1 U P 2 U {z}} , 


Vi€P\j€P*, 
VieP 4 ,j eP\{i}, 
VieP 2 JeP 3 , 

\/i € P 3 ,j € P 3 \{z}. 


Given z E P 4 , j E P 4 , if = 1 is a feasible solution, then to comply with the FIFO 

service order under TSLU operations, the successor of i can only be a pickup vertex in {j}UP 1 UP 2 
or a delivery vertex of a request in P 1 U (P 4 \{j})- Meanwhile, a valid predecessor of j can only be 
a pickup vertex in {z} U P 1 or a delivery vertex of a request in P\(P 3 U {z, j}). Similar reasoning 
can be applied to determine the possible successors of z and predecessors of j for z belonging to 
P 1 and for satisfying CLO service order under TSLU operations. In summary, if arc (z + n, j + n ) 
is included in the routing path, then the set of possible successors of i are obtained as 


Si(i + n,j + n) = {j} U (P 1 ^}) UP 2 Ll{Hn: k E P 2 U P 3 U {z}} , 
Si(i + n, j + n) = {j} U P 1 U P 2 U {k + rz : k E P 1 U (P 4 \{j})} , 
Si(i + n,j + n) = {j} U P 1 U (P 2 \{z}) U {A: + n : k E P 1 U P 4 U {z}} , 
5j(z + n, j + n) = {j} U P 1 U P 2 U {fc + n : k € P 2 U (P 3 \{j})} , 


Vi E P 4 ,j E P 4 , 
VzEP 4 ,jEP 4 \{z}, 
Vz E P 2 ,j € P 3 , 
VzEP 3 ,jEP 3 \{z}, 


and meanwhile the set of possible predecessors of j are 


Vj(i. + n,j + n) = {0} U P 1 U (P 4 \{j}) U {k + n: k E P\{z, j}} , 
Vj(i + n, j + n) = {z} U P 1 U {k + n : k E P\(P 3 U {z, j})} , 
Vj(i + n, j + n) = {0} UP 2 U (P 3 \{j}) U {A: + n : k € P\{z, j}} , 
Pj(z + n, j + n) = {z} U P 2 U {fc + n : k E P\(P 4 U {z, j})} , 


Vz E P 4 ,J G P 4 , 
VzEP 4 ,j EP 4 \{z}, 
Vz E P 2 ,j E P 3 , 

Vz E P 3 ,j E P 3 \{z}. 


With the above insights, four groups of valid inequalities in analog to those in (30)-(33) can 


be obtained to invalidate infeasible predecessors and successors to a given task: 
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Proposition 4. For each of the vertex pairs ieP'u P 4 , j E P 4 \{i}, or i E P 2 U P 3 , j E P 3 \{z} ; 
f/ie following inequalities hold for the PDP-TSLU: 


%ij ^ y %i-\-n,l — 1? 

(36) 

(i,j) 


%ij ^ y %l,j-\-n — 

(37) 



%i+n,j-\-n H - ^ ^ %il ^ 1? 

(38) 

l£Si(i+nj+n) 


%i-\-n,j-\-n + ^ ^ %lj ^ 1* 

(39) 

l£Vj(i+n,j+n) 



The next inequality is adapted from the one proposed for TSPPDF in [22]. It comes from the 
fact that whenever an arc (i, i + n) for any i E P 3 U P 4 is traversed, it follows from the FIFO 
service restriction that the vehicle must arrive empty at i and leaves empty from i + n. 

Proposition 5. For each of the vertex pairs i E P 4 , j E P 4 \{i}, or i E P 3 , j E P 3 \{i}, the 
following inequality holds for the PDP-TSLU: 

Xji + T Xi-\- n j'+n S: lj (40) 

where the variables are zero if their corresponding arcs are not defined in A. 

The following inequality is again adapted from [22], which merges the two inequalities Xij + 
Xj+n.i+n + Xj +7lt i < 1 and x^ + Xj +n j +n + Xj+nj < 1 , for each pair i and j in proper request sets. 

Proposition 6. For each of the vertex pairs i E P 1 U P 4 , j E P 4 \{i} ; or i € P 2 U P 3 , j E P 3 \{i}, 
the following inequality holds for the PDP-TSLU: 


(41) 


where the variables are zero if the corresponding arcs are not defined in A. 


The valid inequalities (26)-(41) will be used to enhance the branch-and-cut algorithm embed¬ 
ded in CPLEX for solving the PDP-TSLU, and their usefulness will be evaluated via numerical 
experiments. 


Remark 3. The above valid inequalities are not exclusive to the PDP-TSLU and other derivations 
are possible. For example, valid inequalities can also be derived from the CFI service restrictions in 
a similar way. Our experiments have shown that they do not contribute much to the computational 
efficiency, and hence are ignored in our studies. For another example, valid inequalities can be 
derived from the geometric restriction of a linear track, known as the “edge degree balance” property 
m ■ The property says that, for an RGV (which starts and ends at the same position) traversing 
any track segment, the number of times it moves towards the right must be equal to that towards 
the left. The property remains true for the PDP-TSLU if a virtual arc is introduced to connect 
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Figure 11: Initial delivery requests and their completion as virtual PD requests. 


the end vertex with the start vertex. Simulations have shown that such valid inequalities hardly 
accelerate the solution process of PDP-TSLJJ and hence are ignored also. 


4. Rolling-horizon approach for handling dynamic PDP-TSLU 

In practice PD requests arrive stochastically and the PDP-TSLU becomes dynamic. Ideally the 
dynamic PDP-TSLU should be solved to optimality at each decision point as a static alternative 
which sequences all PD tasks up to the decision point. In reality, however, a computer has limited 
computational capacity and a suboptimal solution which considers a limited number of PD requests 
has to be explored at each decision point. This motivates us to adopt a rolling-horizon approach 
for handling the dynamic PDP-TSLU. To that end, we need to treat PDP-TSLU with nonzero 
initial load first. 


4-1. Handling nonzero initial conditions 

At the time of recomputing a routing solution, the RGV may contain containers which have 
not been delivered, i.e., ivq > 1. In this case, the re-optimization is subject to a nonzero initial 
condition. The initial load on the RGV may consist of containers that are destined to stations on 


the same or different sides of the track, as depicted in Figure 11, where the impossible types of 
initial requests are avoided by the previous routing solution. 

To recompute a routing decision, extra constraints are required for handling new requests 
without causing conflict with the existing load. Characterizing these constraints and incorporating 
them into the present routing model are thus necessary. The renewed problem can be reformulated 
as a standard one with zero initial load and partially determined arcs, as explained below. 

Let there be no delivery requests at the time when a renewal of routing is triggered, of which 
no,i are destined to north stations and no ,2 (equal to no — 77,0,1 ) to south stations. In the graph 
model, it is feasible to represent the no,i transport requests by no,i virtual PD requests of Type-1, 
each associated with a travel distance equal to the current distance of the RGV to the destined 
station (because visiting a virtual pickup vertex does not introduce a routing cost). Similarly, the 
rest no ,2 transport requests can be represented by no ,2 virtual PD requests of Type-2. See Figure 


11 for a depiction of the virtual PD requests, as indicated by arrows comprised of dashed and solid 


lines. In the meantime, the current position of the RGV is represented by a virtual start vertex 
as directed to one of the virtual pickup vertices with zero routing cost. The order of visiting the 
virtual start and pickup vertices is specified in such a way that yields the initial load setting. 
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The routing model is then augmented to include the initial delivery vertices, the virtual pickup 
and start vertices, and the arcs associated. Denote the set of virtual pickup vertices of Type-1 as 
Pq 1 = {1, 2,..., ?ro,i}, and those of Type-2 as Pq = {no,i + l, no,i + 2,..., no}, and let Po — Pq UPq , 
as a subset of the augmented pickup vertex set P. Without loss of generality, we let the RGV 
start at the virtual start vertex 0 and visit the virtual pickup vertices in the order of 1, 2,, no, 
resulting in the initial load setting. The arcs associated with these virtual vertices ({0} U Po) are 
thus known a priori and the related arc set reduces to 

A 0 = {(k, k + 1) : k = 0, 1,. .., n 0 - 1} 

U {(n 0 , j) : j G P 1 U P 2 U {n 0 ,i + h,n 0 ,2 + h}} 

U {(n 0 , j) : j € P 3 , if n 0 ,i = 0} 

U {(n 0 , j) : j G P 4 , if n 0)2 = 0} , 

where h is a given size of the rolling horizon, satisfying h > hq. All arcs defined in A directing 
to {0} U Po become null since all these vertices have already been visited. For each arc (i,j) in 
the first subset of Ao above, the arc distance is equal to zero; and for each arc (i,j) in the rest 
subsets of Ao above, rij is equal to the distance of the current position of the RGV to vertex j. 

Consequently, the RGV routing problem with nonzero initial load is transformed into a PDP- 
TSLU with partially determined arcs, in which an empty RGV starts at a virtual vertex 0, visits 
the virtual pickup vertices in a predefined order, and then the rest pickup and delivery vertices 
within the rolling horizon by following an order determined by solving the PDP-TSLU. 

4-2. The rolling-horizon approach 

As PD requests arrive stochastically in a complex manner, we adopt a rolling-horizon approach 
to handling the dynamic PDP-TSLU. To make the approach work in real time, we need to address 
two issues: One is to select an appropriate rolling-horizon size and the other is about online 
implementation of the routing decisions. In principle, the rolling horizon should be as long as 
possible in order to avoid myopic decisions, but in the meantime it must be short enough to afford 
real-time decisions. The size of the horizon is thus limited by the computational time required for 
solving a static PDP-TSLU, and the specific size needs to be determined via simulations. Regarding 
the implementation of rolling decisions, two situations have to be handled appropriately: One is 
when re-optimization is demanded while the previous optimization is still in progress; and the other 
is when a new decision is available while the RGV is on the way of executing the previous decision. 
With these situations in mind, the rolling-horizon approach is designed to work as follows. 

Rolling-horizon approach: Whenever a new PD request is issued or there are un-sequenced PD 
tasks in the system, the approach recomputes a solution of sequencing of PD requests by solving 
the (augmented) PDP-TSLU with a given horizon size (if there were so many requests), say, h. 
If there are more than h PD requests to sequence, h of them are selected based on the “Earliest 
Due First” rule while satisfying the FIFO service order at each station. In this way, the RGV 
keeps serving PD tasks as sequenced by the previous decision until a new decision comes out. 
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In addition, executing the first task of each decision is made compulsory and a new decision is 
adopted only after the RGV completes a current task of the previous decision. 

Following this approach, an RGV is able to handle PD requests continuously without obvious 
disruption although the decision is renewed over time. This is true as long as the computation of 
a new decision can be finished before the RGV completes all requests in the previous horizon. 


5. Computational studies 


A toy example is first given to illustrate the modeling and solution of a static PDP-TSLU. 
Then computational results of static PDP-TSLUs with randomly generated requests of various 
sizes are presented to evaluate the usefulness of the valid inequalities for solving the PDP-TSLUs. 
Based on the results, an appropriate horizon size is selected for implementing the rolling-horizon 
approach. Then the rolling-horizon approach is applied to handle dynamic PDP-TSLUs of random 
instances, and the results are compared with those obtained with a typical rule-based method. 

The basic scenario settings of the static and dynamic PDP-TSLUs are as follows. In each 
instance, the PD requests are randomly generated within a work area. Specifically, for each static 
instance a queue of PD requests at any station on both sides of the track are generated with a 
size following a uniform distribution within [0,a], where a is a given integer which is different 
for instances of different sizes. For each dynamic instance, the arrival time of a new PD request 
follows a Poisson distribution and its location follows a uniform distribution within the range of 
stations. The distance between two neighboring stations is treated as a unit distance, and a unit 
service time is assumed for each PD request as equally divided between the pickup and the delivery 
operations. In our studies, an RGV is assumed to have a constant mass equal to two units of load. 
The travel distance of an RGV is defined as the total distance it travels before the last delivery, 


and the energy cost of an RGV is defined as per (23). The friction coefficient of the track is set as 
0.05 and the gravitational acceleration as 9.8 m/s 2 . The problems are coded in C+-1- and solved 
via IBM ILOG CPLEX 12.4 021 which runs on a PC with Intel Core2 Duo CPU @ 2.66 GHz and 


2 GB of RAM. 

5.1. A toy example of static PDP - TSL U 


Consider seven PD requests distributed in nine pairs of face-to-face stations as shown in Figure 
[lj whose graph representation is shown in Figure [l2| In the representation, pickup vertices sit in 
FIFO queues while delivery vertices do not follow any precedence order and are represented to sit 
side by side of each other. The RGV starts from location 3. 

Two scenarios are investigated: (a) the RGV starts with empty load, and (b) the RGV starts 
with two units of load on board. The RGV is simulated with different capacities. The numerical 
results of scenario (a) are summarized in part (a) of Table [I] We observe that the energy cost 
drops considerably once the capacity increases from 1 to 2, but it ceases dropping as the RGV’s 
capacity becomes larger. This is because the service restrictions and the weight-dependent energy 
consumption are to prevent full loading and hence limit the benefit of having a larger capacity. 
The results also indicate that there can be multiple solutions for achieving the same energy cost. 
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(a) RGV starting with empty load. (b) RGV starting with two units of load. 

Figure 12: Graph representations of scenarios with empty and nonempty initial load. 
Table 1. Computational results of the toy example. 


Case 

RGV’s 

capacity 

Energy 

cost 

Travel 

distance 

Completion 

time 

Vertex service order 

(a) 

i 

50.47 

38 

63.78 

4-> 1l->7-> 14->3->10—>5—► 12—> 1 
— >8 — >6 — >13 — >2—>9 


2/3/4/6 

42.63 

30 

55.64 

3->10->4-> 1l->7->5-> 14-> 12—s-1 

— >8 —>6— >2—> 13 —>9 


5 

42.63 

30 

55.64 

l->8->6->2->13->9->3->10^4 

->ll->7->5->14->12 

(b) 

2/3/4/5/6 

55.86 

36 

65.58 

1—>4—>10—>2—>13—>11—>8—>3—>17 
—>12—>5—>14—>6—>15—>9—>7—>18—>16 


In comparison, minimizing travel distance results in higher energy costs for all RGV capacities 
tested, whose values are 44.59, 52.43, 52.43, 48.51 and 52.43 for the capacity equal to 2,3,... ,6, 
respectively (while achieving travel distances of 30, 29, 29, 29 and 29, correspondingly). 

In scenario (b), the RGV starts with two units of load on board, one to be delivered to the 
north station at location 6 and the other to the south station at location 9. The other PD requests 
keep the same as in scenario (a). By introducing a virtual start vertex (0) and two virtual pickup 


vertices (1 and 4), the graph model is augmented as in Figure 12(b) By solving the augmented 
PDP-TSLU with zero initial load and partially determined arcs under five different RGV capacities, 
it yields the computational results summarized in part (b) of Table [I] The routing solutions turn 
out to be the same for the five choices of capacity, resulting in a total energy cost of 55.86. In 
comparison, minimizing travel distance gives total energy costs of 56.84, 62.72, 66.64, 66.64 and 
64.68 for the five capacities, respectively (all achieving a travel distance of 36), which are higher 
than the aforementioned minimum counterparts obtained by solving the PDP-TSLUs. 


5.2. Evaluating the valid inequalities and selecting a size for the rolling horizon 


This subsection evaluates the usefulness of the valid inequalities introduced in Section T3 for 
solving the PDP-TSLU and meanwhile determines an appropriate size for the rolling horizon to 
enable online decisions in a dynamic environment. 


The valid inequalities (25)-(41) are added as redundant constraints to the PDP-TSLU. More 


precisely, the valid inequalities (25) with a subtour size of 2 are added, and the special cases of 
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constraints (26)-(29) are added: 


the constraints (26) with S equal to {i,j}, {i,j + to}, and {i,i + n,j}, and the constraints 


(27) with S equal to {* + n,j + to}, {i,j + n}, and {i,i + n, j + ?z}; 


the constraints (28) and (29) with k = 3: Xi +n j + Xji + Xi t i +n + Xj +rLt i +n + 2 Xj^ +n < 2 and 

x i,i+n T *U+n j'+n T %j+n,i "P %ij T ‘^•Ei,j-\-n — 2. 


These simple constraints were once used in solving TSPPDL and TSPPDF 


In the mean¬ 


time, the valid inequalities (30)-(41) unique to PDP-TSLU each has a polynomial cardinality and 


are added as redundant constraints to the PDP-TSLU directly. 


For convenience of investigation, we classify the valid inequalities (25)-(29) inherited from a 


generic PDP as Group 1, and (30)-(35) derived from LIFO and CLO service restrictions as Group 


2 , and (36)-(41) derived from FIFO and CLO service restrictions as Group 3. The usefulness of 


including various combinations of these three groups of valid inequalities for solving PDP-TSLU 
with/without TW constraints is evaluated via 20 random instances. In each instance, the work 
area consists of 10 pairs of face-to-face stations and the TW constraints (if present) are generated as 
delivery deadlines. The deadline for request i is a random variable following a uniform distribution 
within [T; —15, Tf\. where Tj is the completion time of request i for the same instance in the absence 
of TW constraints. 

The average CPU times to solve the instances are summarized in Table [2j We observe that in 
general, use of any of the three groups of valid inequalities is able to reduce the computational time, 
which is however least likely when valid inequalities of Group 1 are applied. This implies that the 
inequalities inherited from a generic PDP is not as useful as those derived from the specific PDP- 
TSLU, and hence the single inclusion of the Group 1 inequalities is not recommended. Instead, a 
combination of the 1-3 or 2-3 groups of valid inequalities is able to reduce the computational time 
in both simulation cases, either with or without TW constraints. 

The simulation results also indicate that using a rolling-horizon size of eight would allow the 
control system to recompute a routing solution within one minute when the RGV has a capacity of 
2 and about three minutes when the RGV has a capacity of 4. This suggests that a rolling-horizon 
size of eight is reasonable for rendering online routing solutions in a dynamic environment. 


5.3. Evaluating the proposed approach for handling dynamic PDP-TSLU 

Random instances with and without TW constraints are generated to evaluate the rolling- 
horizon approach for treating dynamic PDP-TSLUs. In each instance, 50 PD requests are randomly 
generated in a range of 20 pairs of face-to-face stations. With each request is associated an arrival 
time which follows a Poisson distribution with a mean of 0.5 as simulates the stochastic arrival 
practice. 

In the case with TW constraints, each PD request is associated with a delivery deadline, tight 
or loose. Tight deadlines follow a uniform distribution in the range of [50, 80], and loose deadlines 
follow a uniform distribution in the range of [150, 200]. In each pickup queue the request joining 
first is assigned with a closer deadline. Thus the deadline in each pickup queue is sorted in an 
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Table 2. CPU solution time (in second) when different groups (grps.) of valid inequalities were used. 



i p\ 


7 


8 


9 


10 

Average 


Q 

2 

4 

2 

4 

2 

4 

2 

4 


None 

9.5 

17.7 

42.9 

198.6 

95.1 

438.8 

316.2 

1120.2 

279.9 


Grp. 1 

9.2 

16.9 

38.8 

151.5 

111.1 

658.3 

208.6 

2061.7 

407.0 

Without 

Grp. 2 

9.0 

18.0 

26.6 

131.1 

110.0 

352.1 

226.0 

1179.6 

256.6 

TW cons. 

Grp. 3 

6.0 

18.4 

26.2 

134.3 

87.6 

295.6 

227.0 

1000.9 

224.5 


Grps. 2-3 

6.6 

17.0 

31.3 

132.8 

82.4 

237.9 

272.2 

645.7 

178.2 


Grps. 1-3 

7.5 

16.1 

35.7 

144.2 

95.1 

414.1 

283.9 

972.1 

246.1 


None 

9.5 

19.3 

42.3 

185.5 

149.2 

837.6 

387.9 

2691.6 

540.4 


Grp. 1 

10.8 

21.6 

44.7 

200.1 

207.8 

533.0 

530.1 

2301.6 

481.2 


Grp. 2 

00 

oo 

18.9 

39.3 

201.4 

135.9 

596.8 

400.8 

1801.6 

400.4 

With TW 

Grp. 3 

6.5 

19.4 

31.4 

207.5 

110.7 

600.0 

317.7 

1079.5 

296.6 

cons. 

Grps. 2-3 

7.3 

20.0 

34.6 

165.3 

167.8 

790.4 

286.1 

1776.2 

406.0 


Grps. 1-3 

7.1 

18.5 

37.4 

173.9 

146.6 

668.1 

266.9 

1086.9 

300.7 


increasing order from head to the rear. To simulate schedulable practice, one out of three of the 
PD requests are imposed with tight deadlines and the others with loose ones. 

In each instance, the total energy cost is computed by applying each of the two routing ap¬ 
proaches: the rolling-horizon approach and a rule-based approach. The rolling-horizon approach 
was introduced in Section |4j and is implemented with a rolling-horizon size of eight and by includ¬ 
ing valid inequalities of Groups 1-3 into the problem model. In contrast, the rule-based approach 
is heuristic and it works as follows. 

Rule-based approach: In the absence of TW constraints, all PD requests are queued in each 
station by their arrival times. The pickup service follows the priority rule of “Earliest Arrival 
First”, in which the RGV will first load the container with the earliest arrival time and follow up 
with next arrival container without causing unloading conflict with previously loaded containers, 
until the RGV is full. The RGV will then deliver all loaded containers before picking up any 
new ones. If delivery requests are assigned with TW constraints, the priority rule base is shifted 
from arrival time to deadline to ensure compliance of the TW constraints. The popular heuristics 
of “Earliest Due First” is thus adopted, in which the RGV keeps loading containers with closest 
deadlines as long as no service conflict arises until the RGV is full. Similarly, the RGV will not 
pick up new containers until all loaded ones are delivered. The complete priority rule mimics the 
method used in the current AFHS under consideration, and similar rules had been reported in 
the literature mm • This rule-based approach updates the routing decision whenever a new PD 
request is issued, and the simple algorithm is able to recompute a routing decision in real time for 
a problem having up to fifty requests. 

The two approaches are applied to solving random instances of dynamic PDP-TSLUs without 
and with TW constraints. The average simulation results are summarized in Table [3} We observe 
that, compared to the rule-based approach, the rolling-horizon approach is able to reduce energy 
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Table 3. Total energy costs for operating an RGV. 


Energy cost (without TW constraints) Energy cost (with TW constraints) 


Instance 

Rule-based 

approach 

Rolling- 

horizon 

approach 

Rule-based 

approach 

Rolling- 

horizon 

approach 

Q = 2 

Q = 3 

Q = 2 

O 

II 

CO 

Q = 2 

Q = 3 

Q = 2 

Q = 3 

a 

1335.7 

1222.0 

1067.2 

910.4 

1386.8 

1352.4 

1173.6 

1076.0 

b 

1324.0 

1169.7 

1112.9 

1128.6 

1504.3 

1239.7 

1135.8 

1214.3 

c 

571.3 

583.7 

552.7 

477.3 

592.9 

595.9 

569.4 

505.7 

d 

618.4 

555.7 

516.5 

464.8 

651.7 

607.5 

559.6 

471.5 

e 

641.9 

563.5 

567.4 

540.0 

645.9 

571.3 

639.9 

561.3 

Average 

898.3 

818.9 

7 63.3 

704.2 

956.3 

873.4 

815.7 

765.8 

Percentage saving 


15.03% 

14 . 01 % 



14-7% 

12.32% 


cost by up to 15%. However, the percentage saving decreases as the RGV’s capacity increases 
or when TW constraints are imposed on the PD requests. We also note that the rolling-horizon 
approach has an extra advantage of being able to render feasible decisions in the presence of TW 
constraints, in which case the rule-based approach is often hard to offer. 

Another interesting observation is that, when the RGV’s capacity is increased from 2 to 3, 
both approaches achieve lower energy consumption, but the marginal energy savings are much less 


than the naive estimate of 


«/3—"/ 2 
«/ 2 


x 100% « 33.33%. This occurs because the conflict-free service 
constraints and the weight-dependent energy consumption prevent the RGV exploiting its full 
capacity. The insight can be useful to warehouse designers, helping them determine appropriate 
capacities for RGVs which work in a similar environment by evaluating the trade-off between extra 
construction and operational cost it would take and the actual benefit it would bring. 


6. Conclusions 

This work studied an RGV routing problem in an automated freight handling system. The 
problem was formulated as an MILP that aims to minimize energy consumption for an RGV to 
complete pickup and delivery tasks under conflict-avoidance and time window constraints. The 
energy consumption takes account of routing-dependent gross weight and also dynamics of the 
RGV, and the conflict-avoidance constraints guarantee conflict-free service under two-sided load¬ 
ing/unloading operations. Arc reduction and valid inequalities were exploited from the problem 
structure to enhance the MILP solution process. The static problem model and solution approach 
were integrated with a rolling-horizon approach to handle the dynamic routing problem where air 
cargo enters and departs from the system dynamically in time. Numerical experiments suggest 
that the proposed routing strategy is able to route an RGV to transport cargo with an energy cost 
up to 15% lower than a heuristic method implemented in current practice. 

The current routing problem assumes a single RGV serving the system. In practice, however, 
multiple RGVs may be employed which work on a single track or different tracks that handle 
coupled transport tasks, and meanwhile a container may be relayed to its destination via multiple 
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RGVs. The routing problem thus becomes more complex, and collective routing of the RGVs 
and containers will be required for achieving optimal energy efficiency under service quality and 
feasibility constraints. This constitutes interesting and challenging work for the future research. 
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Appendices 


A. Equivalent forms of the CFI and CLO service constraints 

Lemma 1. The following CFI service constraints are equivalent: 

(See 0b Y, W e P 3 , k e p 1 

<=► Y yh+n< Y £ pl > k e p3 > 

and the following CLO service constraints are equivalent: 

(See 0) Y V* J+n = 0, Vj G P\k G P 1 

i : (i, j+n)eA' 

«=>■ Y Vv - Y Vj G P 1 , fc G Pb 

i:(i,j)€A' i: (i, j-\-n)^A' 

Proof. The proofs of the two equivalences follow similar reasoning. For brevity, we only detail the proof for equiva¬ 
lence between the two CFI service constraints. We prove it by contradiction. 

“=P’: It is sufficient to show that, given an arbitrary pair j G P 3 , k G P 1 , it must have Y+n k+n)eA' Vi k+n — 
Yi:(i,k)eA' Vik■ Suppose this is not true, i.e., there exists a pair j G P 3 , k G P 1 such that Yi : (i,k+n)eA' Vi,k+n > 

E i: (i,k)eA'Vik- This implies that Y i: (i,k+ n )eA'V+k+n = 1 and E i: (i,k)eA' vlk = °- which means that request k is 
picked up before service of request j and delivered during service of request j. Equivalently, this means that request 
j is picked up during the service of request k, i.e., JT. (; j)+A' vti = 1- This is contradictory to the given condition 
and hence proves the sufficiency. 

It is proved alike. Given an arbitrary pair j G P 1 , k G P 3 , suppose J2i (i k)+A r Vik = 1- Thus, 

Y.,i-.(i,k)eA'Vik — Yi:(i,k+n)eA' Vi,k+n’ which means that request k can either be picked up and delivered, or be 
picked up but not delivered during the service of request j. This implies that in certain circumstance the request j 
can be delivered but not picked up during the service of request k, i.e., it is possible to have Yi-(i j+ n )eA' Vi,j+n = 
1 > 0 = Yi (i j)eA' Vi,T This is contradictory to the given condition and hence proves the necessity. □ 


B. Big-M formulation of (12) 

The indicator constraint (|12|) can be linearized by using regular big-M formulation as 


bj > bi + Si + tij - r/ij( 1 - Xij), V(i,j) G A, 
wj > Wi + qj - pij( 1 - Xij), V(i, j) G A, 


(42) 


where rjij and pij are large enough constants such that the right hand sides of the inequalities are lower bounds of 
bj and Wj, respectively. Specific lower bounds can be derived by using the constraints (101, (131, (141 and (201, 
resulting in feasible ijij and pij as summarized in Table 0 where l 0 = q 0 = q2 n +i — 0. The detailed derivation is 
given below. 


28 







Table 4. Feasible constant pairs {r/ij, Pa) (V(i,j) e A) 



i G PU {0} 

i G D 

j G P 

j G DU {2n + 1} ( 

( li Si-\-n H - tij , 

J>i ti,i-\-n “h tij Cj—m Q H - Qj ) 

( J-i—n + tij , Q + qi) 

(li—n + tij — €.j—n , Q + qi + qj) 


Valid reformations of the time and loading consistency constraints in (|12j) require that r/ij > bi + Si + tij — bj and 
pij > Wi+qj — Wj for any (i, j) G A. It is sufficient to set rpj and pij as respective upper bounds of 77 ,. = bi+Si+tij — bj 

13|), (H4| and (201, such upper bounds are derived for four 


and p.. = Wi + qj — Wj for a given (i,j). Based on ( jl0| , 
cases of arcs as follows. 

Case a: i G P U {0}, j G P. It follows that 77 .. < b i+n — ti^+ n + tij — bj < /, — Si+„ — ti t i+ n + tij, and that 
P.. <Wi + Pj ~ Qj < Q- 

— L J 

Case b: i G Pu{0}, j G DU{2n+l}. It follows that 77 .. < bi +n — ti,i+ n +tij — bj < h - Si+ n — ti,i+ n +Uj - ej-„, 


and that p.. < Wi + qj < Q + qj. 

—' l J 

Case c: i G D, j G P. It follows that 77.. < U- n + tij — bj < U- n + tij, and that p .. < Wi + qj — qj < Q + qi. 

—tj —tj 

Case d\ i G D, j G D U {2n + 1}. It follows that 77 ,, < U- n + tij — bj < U- n + tij — ej- n , and that 

—'tj 

P-- < Wi + qj < Q + qi + qj. 

—' l 3 
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